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Jf I • Abstract. Quantum interfaces between polarized atomic ensembles and coherent 

states of light, applied recently to manipulate bipartite and multipartite entanglement, 
are revisited by means of a continuous-variable formalism. The explicit use of the 
continuous-variable formalism facilitates significantly the analysis of entanglement 
between different modes, reducing it to the study of the properties of a final covariance 

£h ' matrix which can be found analytically. Furthermore, it allows to study matter-light 

interfaces for mixed states, adapting the formalism to the experimental situations 

qh| in which the initial prepared Gaussian states are, unavoidably, affected by a certain 

amount of noise. A multipartite scenario, leading to the generation of macroscopic 
^ , cluster states is presented and analyzed in detail within this formalism. 
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1. Introduction 

A strongly polarized macroscopic atomic ensemble can be described disregarding its 
^ individual atomic components and using collective spin variables instead. In this 

scenario, the interaction of an atomic polarized ensemble with off-resonant linearly 
polarized light resolves, at a classical level, into a Faraday rotation. The latter refers 
to the rotation experienced by the polarization of light propagating inside a magnetic 
medium (polarized atomic ensemble). Moreover, this interaction can also lead to an 
exchange of quantum fluctuations between matter and light, i.e. a quantum interface. 
Seminal results exploiting such interface are the generation of spin squeezed states [I], 
and the entanglement between two macroscopic spatially-separated atomic ensembles 
[2]. In the first case, squeezing was produced by the interaction of an atomic sample with 
a squeezed state of light, whereas the second task was achieved by propagating a single 
laser bram through both samples. The above protocols demand a final measurement on 
the light (homodyne detection) which projects the state of the atomic ensembles into 
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a desired final state. This mapping of fluctuations between different physical systems 
provides a powerful tool to design quantum correlations. The potential applications 
of such a genuine quantum instrument have just started. Among them the use as a 
spectroscopic method for strongly correlated ultracold atomic systems [3J. Furthermore, 
if the incident light beam is spatially tailored, the interface provides spatial resolution. 
This characteristic should permit to directly observe correlations without the need of 
individual spin addressing jIHH], or even provide a direct measure of the order parameter 
of some exotic quantum states of matter [7]. Conversely, a quantum interface between a 
strongly correlated system and light should allow to map the correlations of the matter 
onto the light, leading to the generation of highly non classical states of light. 

Here we address matter-light interfaces in the framework of continuous-variable 
(CV) formalism which strongly facilitates the manipulation and verification of 
entanglement, both bipartite and multipartite, reducing it to the study of the properties 
of a final covariance matrix. Such formalism is especially well suited for mixed states 
since they are treated on an equal footing as pure states. The formalism of continuous 
variables arises naturally if the initial states of atomic ensembles and light are Gaussian. 
Then, if the interaction between light and matter is bilinear, the Gaussian character of 
the subsystems is preserved. Unitary evolutions of such form correspond to symplectic 
transformations of the covariance matrix of the composite system (matter and light). 
The final projective measurement on the light, homodyne detection, is also a Gaussian 
operation and it is easily performed at the level of the final covariance matrix. Therefore, 
all the necessary tools are at hand. 

The paper is organized as follows. In section 2, we briefly state the problem 
under consideration, starting from the description of the matter and light as 
continuous-variable systems. Then using the Quantum Non Demolition (QND) effective 
Hamiltonian of the matter-light interaction we derive the evolution equations (for a 
detailed derivation of such Hamiltonian the reader is referred to [8]). In section 3, we 
derive explicitly all the steps of the matter-light interface in terms of covariance matrices 
and symplectic transformations [9l[10]. In section 4, we demonstrate the suitability of 
the formalism by analyzing in detail two examples. The first one corresponds to the 
well known EPR bipartite entanglement between two macroscopic atomic ensembles, 
experimentally demonstrated in [2]. We consider a more general scenario in which 
the atomic samples are initially in thermal states with different amount of noise. In 
such a case entanglement between the atomic samples is produced only above a given 
value of the matter-light coupling constant. By exploiting different geometries [JTJ, the 
interaction with a second light beam permits to design optimally the final entanglement 
between the samples and even to delete the entangling action of the first beam. The 
second case we treat, much more involved, analyzes the generation of cluster states 
using matter-light interfaces and demonstrates the full power of this formalism in the 
multipartite scenario. In section 5, we present our conclusions. 
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2. The Faraday Interaction 

The physical system we analyze consists of several atomic ensembles and light beams, 
the latter playing the role of information carriers between the atomic samples. At a time, 
only a single light beam interacts with the atomic ensembles. After the interaction, the 
light beam is measured. 

Each atomic ensemble contains a sufficiently large number, JV at , of noninteracting 
alkali atoms with individual total angular momentum F. The ensemble is described by 
its collective angular momentum J = ( J x , J y , J z ), where J k = Yli=i Fk,i (k = x, y, z). All 
atoms are assumed to be polarized along the x direction, which corresponds to preparing 
them in a particular hyperfine state \F,mp). In such a situation, fluctuations in the J x 
component of the collective spin are very low, allowing for treating this variable as a 
classical number J x ps (J x ) = HJ X = h~N at F. As a consequence, the quantum character 
of the collective spin is preserved in the orthogonal spin components, which have a 
zero mean, but non-zero fluctuations. By appropriate normalization they are made 
to fulfil the canonical commutation relation, J y /y/hJ x , J z j ' \J%J X = ih. To stress the 
continuous variable character of the system, we rename the above variables as "position" 
and "momentum" : 

X A - 



\fh\L x 

\fhT x 

and from now on use only the canonical variables xa,Pa to refer to the atomic sample, 
where the subindex A stands for atomic ensemble. Later on when we deal with few 
ensembles we will use the notation XA,mPA,n to refer to the nth atomic sample. 

On the other hand, the light is taken to be out of resonance from any relevant 
atomic transition and linearly polarized along, say, the x-direction. We use the Stokes 
vector description s = (s x , s y , s z ) of light polarization. The components §k (k = x, y, z) 
correspond to the differences between the number of photons (per unit time) with x and 
y linear polarizations, ±7r/4 linear polarizations and the two circular polarizations, i.e. 
' h h 

h,. , , h 



s„ = -in. 



n\) = o(4oi/ + ojo*)> ( 2 ) 



[n - n ) = —{&la y - a\a 2 



'v 9 2 y 

h,^ , h 

2(n -n ) = -, 

These allow for the microscopic description of the interaction with atoms, however, 
effectively only the following macroscopic observables will be relevant: S^ = j Q §k(t)dt, 
where T is the duration of the light pulse. So defined operators obey standard angular 
momentum commutation rules. The assumption of linear polarization along direction 
x allows for the approximation S x ~ (S x ) = N p hh/2. Once more, the remaining 
orthogonal components S y and S z are appropriately rescaled in order to make them to 
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fulfill the canonical commutation rule, 
correspondence equivalent to equation 

XL ~ 



S y /y/hS x , S z /y/hS x = ih. Straightforwardly, a 
[T]) arises: 



\/hSx 



PL = -7=, (3) 

y/hb x 

which allows to treat the light polarization degrees of freedom on the same footing 
as the atomic variables. Notice, that while only one light beam is the carrier which 
entangles the atomic ensembles, a secondary light beam will be normally needed to 
verify entanglement [8|ITT]. 

In the situation in which a light beam propagates in the YZ plain and passes 
through a single ensemble at angle a with respect to direction z, the atom-light 
interaction can be approximated to the following QND effective Hamiltonian (see [8] 
for a detailed derivation) 

H int ( a ) = ~ ifVL (Pa cos a + x A sin a) . (4) 

The parameter k is the coupling constant with the dimension of the inverse of 
an action. Notice that such Hamiltonian leads to a bilinear coupling between the 
Stokes operator and the collective atomic spin operators. Evolution can be calculated 
through the Heisenberg equation for the atoms and using Maxwell-Bloch equation for 
light, neglecting retardation effect. The variables characterizing the composite system 
transform according to the following equations ( [8] and references therein): 

x°f = ^-Kf, 1 l cos a, (5 a) 

PT = ^ + /^ n sina, (56) 

x° L ut = x l £ - k{Pa cos a + x { % sin a) , (5c) 

PT=fl- M 

The above equations can be straightforwardly generalized to the case in which a single 
light beam (xl,Pl) propagates through many samples shining at the nth sample at a 
certain angle a n . 

Due to the strong polarization constraint the initial states of the atomic ensembles 
as well as the one of light can be treated as Gaussian modes. On the other hand, the 
interaction is a bilinear coupling between the Stokes operator and the collective atomic 
spin operator, thus, it preserves the Gaussian character of the initial modes and can be 
interpreted as a Gaussian interaction between two bosonic modes. These facts enable 
us to tackle the quantum interface within a CV formalism. 

3. The matter-light interface in the CV formalism 

We start by reviewing the most basic concepts needed to describe Gaussian continuous- 
variable systems. For further reading, the reader is referred to [i~2HT4"] and references 
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therein. For a general quantum system of iV pairs of canonical degrees of freedom 
("position" and "momentum"), the commutation relations fulfilled by the canonical 
coordinates R = (xi,pi,...,Xn,Pn) can be represented in a matrix form by the 
symplectic matrix J N : [R\, Rj] = ift(i7iv)y> hj — 1> • • • > 2iV, where 

J N = 0J, J=f° i jY (6) 

Gaussian states are, by definition, fully described by the first and second moments of the 
canonical coordinates. Hence, rather than describing them by their infinite-dimensional 
density matrix g, one can use the Wigner function representation 

" f W = ?f^»PHC-<0V , (C-<0], P) 

which is a function of the first moments through the displacement vector d, and of the 
second moments through the covariance matrix 7, defined as: 

di = Tr(gRi), 7^ = Tr(g{Ri - d i} Rj - dj}). (8) 

The variable ( = (xi,pi, . . . ,Xn,Pn) is a real phase space vector with probability 
distribution given by the Wigner function. The covariance matrix corresponding to 
a quantum state must fulfill the positivity condition 

7 + i Jn > 0. (9) 

In the particular case of a physical system consisting of several atomic ensembles and 
single light beam the most general covariance matrix takes the form 

( 1 A c 
1 ~\C T 7 L 

where the submatrix 7 L corresponding to light mode, j A the atomic ensembles, and C 
accounts for the matter light correlations. 

If a Gaussian state undergoes a unitary evolution preserving its Gaussian character, 
which is the case in the physical systems under consideration, then the respective 
transformation at the level of the covariance matrix is represented by a symplectic 
matrix S acting as 

7out = S T j in S. (11) 

We illustrate how to reconstruct the evolution of the covariance matrix from 
the propagation equations (l5aj) - (|5dj) . One notes that the variables describing the 
system after interaction are expressed as a linear combination of the initial ones. 
Let us denote by K the following linear transformation K : (x 



(10) 



out 
A,ni 



out 



Pa 



~ out "out \ T 
X L iPl I 



K ( x a m Pa m z™ i P 1 l) T i which can be straightforwardly obtained from equations 



and for a single atomic mode reads: 
/ x°f \ ( 

"out 

Pa 



x 



out 



V PT ) 



1 








— k cos a ' 





1 





Ksina 


Ksina 


— ncosa 


1 














1 I 



( $1 \ 

P>A 
"in 
X L 

\Pl j 



(12) 
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Since the interaction Hamiltonian is bilinear, the matrix K can be directly applied to a 
phase space vector ( and correspondingly to the covariance matrix, however the sign of 
the coupling constant k should be changed. This is because the phase space variables 
undergo the Schrodinger evolution, whereas the quadratures transform according to the 
Heisenberg picture. Therefore, we define K = K\ K _^, K ^ which we apply to the phase 
space vector and covariance matrix as 

'soutTin S>out = Cin-^ 7in "■ Cin 

= Cin7outCin) (13 J 

leading to S = (K 7 )^ 1 . The above formalism has been explicitly developed for a single 
sample and a single beam, but it easily generalizes to an arbitrary number of atomic 
samples and beams, as well as to various geometrical settings. 

Finally, the last ingredient essential to describe the matter-light interface at the 
level of the covariance matrix is the homodyne detection of light. Assuming a zero 
initial displacement and covariance matrix of the form ( TTOj) . the measurement of the 
quadrature xl with outcome xl leaves the atomic system in a state described by a 
covariance matrix 



/' = 1 A -C(X 1 L X)- 1 C T , (14) 

and displacement 

d A = C(X 1 L X)- l (S: L ,0), (15) 

where the inverse is understood as an inverse on the support whenever the matrix is not 
of full rank and X is a diagonal matrix with the same dimension as 7 L with diagonal 
entries (1, 0, 1, 0, . . . , 1, 0). 

Inherent to the matter-light interface is the analysis of the remaining atomic 
samples once the light beam has been measured , in order to check the correlations 
induced between them. With the CV formalism we have access to the complete 
information provided by the covariance matrix and displacement vector after interaction. 
This makes verification of entanglement amenable to covariance matrix criteria (see 
also 0EQ]). 

A structural separability test, which can be only applied when a full covariance 
matrix is available, is the positive partial transposition (PPT) test [T6l[TT] . For 
continuous-variable systems, it corresponds to partial time reversal of the covariance 
matrix [18], i.e. the change of the sign of the momentum for chosen modes. If 
the partially time reversed covariance matrix does not fulfill the positivity condition 
(J9]), the corresponding state is entangled. This test, however, checks only for the 
bipartite entanglement. For Gaussian states this criterion is necessary and sufficient 
for entanglement of 1 x N modes. 

Experimentally more convenient separability test is based on variances of collective 
observables. It was provided for two-mode states in [Tj5] and generalized to many-mode 
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states in [20]. It states that if an N mode state is separable, then the sum of the 
variances of the following operators: 

u = h\X\ + . . . + h^XN 

v = gipi + • • • + 9nPn (16) 

is bounded from below by a function of the coefficients h\, . . . ,h]sr,gi, ■ ■ ■ ,9n- 
Mathematically, the inequality is expressed as 

(An) 2 + (At;) 2 > f(h h . . . , h N ,g h . . . , g N )h, (17) 



where 



f(h l ,...,h N ,g 1 ,...,g N ) 



hgi + ^ h r9r 



+ 



h m g m + 2j ^s9i 



ser 



(18) 



r£l 

In the above formula the two modes, I and m, are distinguished and the remaining ones 
are grouped into two disjoined sets I and I'. The criterion (j!7p holds for all bipartite 
splittings of a state defined by the sets of indices {1} U I and {m} U I'. For two mode 
states, the criterion becomes a necessary and sufficient entanglement test, however only 
after the state is transformed into its standard form by local operations [19J. This local 
transformations, however, are determined by the form of the covariance matrix. In this 
sense, the knowledge of the full covariance matrix is essential in order to determine 
whether the state is entangled. Since in experiment we usually do not have access to 
this information, we cannot assume that the criterion decides unambiguously about 
separability. 

4. Symplectic and Covariance matrices for bipartite and multipartite 
entanglement 

Let us illustrate the versatility of the formalism in two distinct examples, the first 
one deals with entanglement of two atomic samples generated and manipulated by 
interaction and measurement of a light pulse. In the second one, with larger number 
of atomic samples, we want to show how a cluster state can be generated using the 
atom-light interface. 

4.1. Bipartite entanglement for mixed states: Generation and verification 

In the seminal work of Polzik and coworkers [2], the entanglement between two spatially 
separated atomic samples was generated in the experimental setup schematically shown 
in figure [TJ In such setup, both light and atomic samples were strongly polarized along 
the x-direction while light propagated along the z-direction (a = 0). The atomic 
ensembles were previously addressed with local magnetic fields oriented in opposite 
directions, enabling to use a single light beam for generation of EPR entanglement 
and another one for verification, however, between new time-integrated variables. It 
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-B B 



z 




Figure 1. Schematically presented experimental setup used in [2] to generate 
and verify the presence of entanglement between two macroscopic atomic samples. 
The external magnetic field makes feasible the measurement of the two transverse 
components of the spin with a single light passage. 



is straightforward to generalize the equations of motion (|5al 
interacting with light according to the Hamiltonian 



H 



(2) 



K 

T 



Pl(pa,i +Pa,2) 



5dj) to two samples 



(19) 



For example, in this case the canonical "momentum" variable for atoms, pA, n (n = 1, 2), 
is preserved [compare with equation (J5jjj) ], while the canonical "position" of light after 
the interaction carries information about the sum of atomic momenta 



x 



out 



X, 



. h -K (p A ,l + PAfl) ■ (20) 

Entanglement between the atomic samples is established as soon as the xf^ component 
of light is measured. Moreover, it should be emphasized that entanglement is generated 
independently of the outcome of the measurement, nevertheless the measurement result 
indicates the displacement of the state. 

The setup described above needs some modifications if the individual magnetic field 
addressing is impossible. The way to overcome this problem was shown in [TTj, and is 
summarized in figures |2i and [2b. 

Here we give a detailed description of this setup at the level of covariance matrix. 
Moreover, we further assume that the initial state of the two samples is not a minimum- 
fluctuation vacuum state as in [8], but a general thermal state. Under such assumptions, 
the initial state of the composite system is given by the following covariance matrix for 
atoms and light 7 in = nilLf © n 2 l^ © if, where the identity 1 2 stands for a single 
mode and parameters ni,n2 are related to temperature by rij = 1/ ta.nh[hw/ (2k sTi)] 
{i = 1,2). The symplectic matrix, S^t, describing the interaction of light passing 
through the samples at zero angle (figure [2ji) is given by 



S 



int 



/ 1 


\ 


10 


K 


10 





1 


K 





1 


\ K K 


1/ 



(21) 
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a) 



®X 



■^z 




Figure 2. The sketch of the setups using the geometrical approach to generate 
and manipulate bipartite entanglement. The interaction between the light beam and 
atomic samples followed by the measurement introduces squeezing in a) pa,i +Pa,2 b) 
xaa - xa.2 and c) Xaa + xa.2- 



thus, the covariance matrix after the interaction takes the form expressed in equation 
(HTD 



7out 



/ m + k 2 





K 2 








k\ 





m 








n\K 





K 2 





n 2 + k 2 








K 











n 2 


n 2 K, 








n\K 





n 2 n 


1 + nxK 2 + n 2 K 2 





\ K 





K 








1/ 



(22) 



Both modes, representing the samples, are entangled with light, however their reduced 
state is separable as one can check applying PPT criterion to the covariance matrix 
of the upper-left block matrix. Entanglement between atomic samples is not produced 
until one measures a quadrature of light. To demonstrate this, we apply a homodyne 
measurement of the light mode [see (TH1) and (TT5|) ]. Assuming the measurement outcome 
Xla obtaining the covariance matrix describing the final state of the samples 



7fin 



/ n x + k 2 


2 



V 



K 







nira2K 2 +rti 
(ni+n 2 )K 2 +l 



n±n2K 2 
(ni+n 2 )K 2 +l 





n 2 + k 2 









\ 



nxn^K, 



(ni+n 2 )K 2 + l 



niri2K 2 +n2 
(ni + n 2 )K 2 +l 



(23) 



and the displacement of the final state is 

Notice that the covariance matrix is independent of the measurement outcome, but the 
latter is present in the displacement vector. 
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In order to verify that the final state of atomic samples is entangled we use the 
separability criterion based on the variances of the two commuting variables 



A(\MPa,1 + jPA,2) 



-i 2 



+ 



A(\X\x A ,l - T^Afi) 



1 2 



> in. (25) 



We emphasize, however, that due to the fact that in the analyzed setup the measurement 
possibilities are limited, the only experimentally feasible case is |A| = 1. The way to 
measure such combination of variances was described in detail in [TT] , therefore we will 
not recall it here. Also the final state is not, in general, in its standard form, therefore, 
the measurement gives only a sufficient condition for separability. 

In the following analysis we restrict to a single inequality involving the collective 
observables which are squeezed during entanglement generation, i.e. the one for A = 1. 
From the final covariance matrix (1231) one can directly compute the variances of the 
collective atomic spin: 

1 2 1 

- [A(pA,l + PAfl)} = ~ (7fin,22 + 7fin,44 + ^^24) 

m + n 2 



2(n 1 + n 2 )/t 2 + 2' 



1 2 1 

- [A{XA,1 - XA,2)\ = 2(7fin,H + 7fin,33 ~ 27fi n ,i 3 ) 

= \{n 1 + n 2 ). (26) 

Substituting the obtained variances into the separability criterion we obtain 

rw- , . v,2_, r/ w~ * m2 (gi + n 2 ) (wi + wg) /t 2 + 2 

A (p A ,l +PA,2)\ + A (x A ,i - x A ,2)\ = 7 — - — v 2 - , 27 

2 [m + n 2 ) k 2 + 1 

which violates the bound 2h~ for specific values of m, n 2 , and k, i.e. for those that fulfill 

y^y/nx + n 2 - 2 
k > . (28) 

V (-ni - n 2 + 4) (m + ri 2 ) 

The criterion (T25|) for A = 1 does not detect all entanglement. We compare the set that 
is not detected by the inequality to the set of separable states found using the PPT 
criterion. The respective ranges of parameters K,ni,n 2 are depicted in figures E^, and 

Eb. 

The comparison shows that to generate entanglement between thermal states, much 
stronger coupling is necessary. Also, the simple criterion based on the variances cannot 
detect entanglement for n\ + n 2 > 4, even though it is clearly present in the system. 

The procedure summarized in figure |2b introduces more entanglement in the system 
squeezing the collective variable x A ,i — %a,2- The computation is analogous, therefore 
will not be repeated. The resulting state is characterized by the following variances: 

s [A( ^ + p " )|2 = s |A(i '« " ^ f = 2 (ni "'^ + 2 - (29) 



The condition for the inequality (1251) for A = 1 to be violated is 



K y '/ 2n, + 2„, ■ (30) 
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Figure 3. Comparison of the sets of parameters k, 711,712 for which the states are 
entangled [figures a) and c)], with the ranges of parameters for which they are detected 
by the spin variance inequality ([25)1 with A = 1 [figures b) and d)] . In figures a) and b) 
we consider the state produced in the setup from figure [2k-, whereas in figures c) and 
d) the one produced in two steps schematically depicted in figures [5^ andj^h- 



The set of parameters k, nl, nl for which the final state is detected by the inequality is 
depicted in figure [3b. For comparison the set of separable (equivalently PPT) is shown 
in figure EH. Again the inequality does not detect all entangled states, however is more 
efficient than in the first case. 

Interestingly enough, our geometrical approach, making the light impinging on 
each atomic sample at a given angle a,, also opens the possibility of deleting all the 
entanglement created by the first light beam, by interaction with a second light beam 
of an appropriate intensity. Notice that the entanglement procedure is intrinsically 
irreversible because it involves a projective measurement, so coming "deterministically" 
back to the initial state is not obvious. In [2T1I22], a quantum erasing scheme in 
continuous-variable systems was proposed. The measurement of the meter coordinate 
entangled with the quantum system causes a back-action. The authors shown that 
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it is possible to erase the action of the measurement and restore the original state of 
the system. Here we are interested in deleting the measurement induced entanglement 
between two atomic samples, exploiting the squeezing and antisqueezing effects produced 
by the laser beams. Again, a CV formalism greatly simplifies the analysis and indicates 
in which way the atomic quadratures, squeezed and anti-squeezed by the first light 
beam, can be restored with the help of a second light pulse. 

We begin with the state produced in the first part of the section, represented by the 
covariance matrix ( |23|) . We will perform this part of the analysis only for the initially 
vacuum state, i.e. n\ = n 2 = 1. Evolution due to the interaction with the light beam 
impinging on each atomic sample at ct\ = a 2 = 7r/2, schematically represented in figure 
[2b is reproduced by the following symplectic matrix acting on 7 fin © 1 2 : 



Ccrascr 
D int 



(1 








V o \ 





1 











1 





7] 








1 














1 


\o 


-7] 


-n 


l) 



(31) 



In general, the second light beam may have different properties, hence we use a different 
coupling constant n. After the interaction and measurement of light the state of atomic 
ensembles is characterized by the following covariance matrix: 



/ rj 2 (2n?+l)+K 2 +l 



72nd 



2^2(2 


K 2 + l) + l 







K 2 —ri* 


(2« 2 +l) 







K 2 -r) 2 {2K 2 +l) 



\ 



2?7 2 (2k 2 +1)+1 





2t? 2 (2k 2 +1)+1 
rj 2 (2k 2 +\)+k 2 +1 

2k 2 +1 

»? 2 (2k 2 +1)+k 2 +1 







k 2 -»; 2 (2k 2 +i) 
2k 2 +1 



2ti 2 (2k 2 +1) + 1 







K 2 -r7 2 (2K 2 + l) 
2k 2 + 1 



r ? 2 (2K 2 +l)+K 2 +l 
2k 2 + 1 



\ 



/ 



For the specifically adjusted values of n and n, i.e. rf = k 2 /(1 + 2/t 2 ), the atomic 
ensembles come back to the initial coherent state, at least at the level of variances. 
Let us have a look at the displacement vector. Assuming that the outcome of the 
measurement leading to the generation of entangled state ( |23l) is xl,i (see equation 
( I!M|) ). and the outcome of the measurement erasing the correlations is x^, the total 
displacement is given by 



d, 



K c 



2k 2 + 1 



X L,2, 



2k 2 + 1 



X L,1, 



K 



-X L ,2 



2k 2 + 1 ' 2k 2 + 1 



X L A 



(33) 



Hence the final state is equivalent to the initial one only up to the displacement vector. 



.(32) 



4-2. Multipartite entanglement: Generation of cluster states using atom-light 
interfaces. 

In [23] a class of A-qubit quantum states generated in two dimensional arrays of qubits 
using an Ising-type interaction was presented; these are the so-called cluster states. 
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Figure 4. Schematically depicted setups leading to generation of the four-partite 
linear cluster state. Figure a) shows the graph representing the state structure. In 
figures b)-e) the sequence of beams introducing the CV cluster-like corelations between 
atomic samples are depicted: b) p' A 1 — x' A 2 , c) p' A 2 — x' A 1 — x' A 3 , d) p' A 3 — x' A 2 — x' A 4 
and e) #4 4 -£^,3. 



Using their scalability properties, Briegel et al. proposed a scheme for a one-way 
quantum computation [23J|25]. There, two-dimensional cluster states are the entire 
resource, while computation consists of a sequence of local projective measurements. 
Furthermore, any two-dimensional cluster state provides a universal quantum computer 
since it has been proved that any unitary quantum logic network can be efficiently 
simulated within this scheme. In this sense cluster states, can be regarded as a resource 
able to generate any type of multi-qubit entanglement by means of two body interactions. 
We present here, how the continuous- variable cluster-like states [251 EZ] can be 
generated within the analyzed atom-light interface. We associate the modes of the N s - 
mode system with the vertices of a graph G. The edges between the vertices define 
the notion of nearest neighbourhood. By iV a we denote the set of nearest neighbours 
of vertex a. A cluster is a connected graph. For angular momentum variables, cluster 
states are defined only asymptotically as those with infinite squeezing in the variables 

Pa~Yl ^ ( 34 ) 

beN a 

for all a G G. We talk about the clusterlike states when squeezing is finite. 

In [UJ we have shown how to generate a simple cluster-like state between four 
samples (figure Hk), however in a rotated frame. For each mode we introduce new 
variables: 

%A,n = ~7^(^A,n —pA,n), V A,n = ~7^{ x A,n +pA,n)- (35) 

The setup is summarized in figure HI 
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The state is generated in four, commuting steps. We will denote an initial state of 
light and atoms at each step by 7^ , a state after interaction by 7^., and the state after 



U) 



the measurement of light by 7^ j 



4. Note that 7 



(i+i) 



7^ © 1 2 . 



,(i) 



We assume that initially both the atomic samples and light are in a vacuum state 
7 in ' = I-8+2 • The symplectic matrix (in a rotated frame) reproducing the effect of 
interaction depicted in figure Hb is 

\ 



s\ 



(1) 



( l ° 








1 


K 






-K 












1 








00 

k ° 


1 


^ K 


1 



(36) 



/ 

whereas the interaction from figure Hfc is reproduced with symplectic matrix 

\ 



S\' 



(2) 



f l ° 





-K 


1 
















K 





1 


-K 









000000 


1 


V K 


K K 


1 



(37) 



/ 

The other two symplectic matrices can be easily written down basing on Sj and Sj . 
The final state after the sequence of four atom-light interactions followed by the 
measurement of light is characterized by the covariance matrix 

/ A C E F \ 
C B D E 



(4) 
7fin 



E 
F 



D 

E 



B 

C 



C 
A 



(3* 



where A, B, C, D, E and F are the following 2x2 matrices: 

5k 6 +8k 4 +5k 2 +1 
1 



A 



B 



C 



5k 4 +5k 2 +1 



5k 6 +7k 4 +4k 2 + 1 
5k 4 +5k 2 +1 




k 2 (5k 4 +7k 2 +2) 
5k 4 +5k 2 +1 





5k 6 +7k 4 +5k 2 +I 



5k 4 +5k 2 +1 



10k 6 + 13k 4 +6k 2 + 1 
5k 4 +5k 2 +1 

k 2 (5k 4 +7k 2 +2) 

5k 4 +5k 2 +1 





D 



E 







k 2 (5k 4 +6k 2 +2) 

5k 4 +5k 2 +1 
^.4 Jl 



5k 4 +5k 2 +1 





k 2 (5k 4 +6k 2 +2) 
5k 4 +5k 2 +1 











5k 6 +6k 4 +k 2 
5k 4 +5k 2 +1 

!-4 



(39) 



5k 4 +5k 2 +1 



5k 4 +5k 2 +1 
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In order to check the separability of the produced state, we apply the PPT criterion. 
The partial time reversal of covariance matrix ( 158)) with respect to all cuts is negative. 
Therefore, the state is fully inseparable for all values of k. 

Applying a less effective, but experimentally convenient separability test of the 
form (ITT)) one detects full inseparability above some threshold value of k. The three 
inequalities which constitute a necessary separability criterion for this particular cluster 
state contain variances of the variables squeezed by the interaction and read 



A 1 = [A(p[ - x' 2 )f + [A(p> 2 - x\ - x' 3 )f > 2h (40a) 

A 2 = [A(& - x' 2 - x' A )f + [A(p' 2 - x\ - x' 3 )} 2 > 2h (406) 

A 3 = [A(p' 2 - x[ - x' 3 )f + [A(p' 4 - x' 3 )f > 2h. (40c) 

We determine the variances explicitly from the covariance matrix (138)) : 



\[A{p'^-x^)] 2 = l -A 22+ l - Bll -C 2l 

5k 2 + 2 



10/t 4 + 10k 2 + 2' 



(41a) 



l - [A(/ 2 ) - x'« - x'^)] 2 = \b 22 + l -A xl + l -B u - C 12 - D 21 + E u 



5k 



2 



10k 4 + 10k 2 + 2 : 



(416) 



\ [A(/ 3 ) - x'^ - x'^)] 2 = \b 22 + l ~B n + \A X1 - D 12 - C 21 + E u 



h i v ,j 2 2 2 

5/« 2 + 3 



10k 4 + 10k 2 - 

l -[A{p'^-x'^)] 2 =\A 22 + \Bn-C 12 

5k 2 + 2 



(41c) 



, -• (41 d) 

10k 4 + 10k 2 + 2 v ; 

The dependence of the left-hand sides of inequalities (140 aty on the coupling constant is 
depicted in figure |5l The inequalities (140 aft and (140 c\i are violated for k > 0.3 and the 
inequality (140 6j) is violated for k > 0.4. 



5. Summary 

We have addressed the problem of manipulation of entanglement between atomic 
ensembles using Faraday quantum atom-light interface with a continuous-variable 
formalism. Both the atomic ensembles and light can be prepared in initial Gaussian 
state and remain Gaussian after all the steps leading to generation of entanglement. 
This allows us to describe the whole process using the covariance matrix. After a 
general introduction of the mentods we have applied it to derive bounds on the strength 
of the interaction leading to entanglement between atomic ensembles prepared initially 
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H 



Figure 5. Dashed line - equation (141 q[) . dashed-dotted line - equation (|41b[) . purple 
equation (|406[) . orange - equation (|40al) and (|40cp 



in mixed states. Then we have shown how the CV formalism facilitates the analysis of 
more complex systems. In particular we have addressed the problem of generation of 
the CV cluster states. 
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